Effective field theory as the bridge between 
lattice QCD and nuclear physics 
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Abstract. A confluence of theoretical and technological developments are beginning to make 
possible contributions to nuclear physics from lattice QCD. Effective field theory plays a critical 
role in these advances. I give several examples. 
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INTRODUCTION 

While it is unrealistic to expect to see a solution of the structure of a uranium nucleus 
from QCD within our lifetimes, it is not unreasonable to predict that lattice QCD 
will make significant contributions to nuclear physics over the next couple of decades. 
Simultaneous progress in computer technology, computational algorithms, and advances 
in theory have made it feasible to begin such a program in lattice nuclear physics now. 

The limitation one faces is the computational cost of a realistic simulation. L. 
Giusti presented the following formula at Lattice '06 for the cost (in Tflops- 
yrs) for generating gauge field configurations with dynamical Wilson fermions 
( |http://www.physics.arizona.edu/lattice06/| ): 



Cost -0.15 



#configs 
1000 



m 



L20MeVJ 



V 



32 fm^ 



0.08 fm 



(1) 



Here a is the lattice spacing, niq is the light quark mass, V is the lattice volume. 

Significant advances in algorithms have occurred in recent years, with the discovery 
in the 1990's of how to simulate chiral fermions fl','2!], and the improvement of methods 
for including light dynamical fermions (for example, the power of the mass dependence 
in the above formula has dropped from ^ to m^^ since the development of algorithms 
in refs. [3, 4]). Technological advances continue unabated, and machines currently exist 
operating in the 10^ Tflops range, and Pfiops computing will exist before long. 

Nevertheless, technology plus algorithms do not by themselves add up to advances 
in nuclear theory in the near future because of the daunting computation costs of a 
realistic simulation. To avoid the disadvantages of non-chiral lattice fermions, such as 
the Wilson formulation, one should use domain wall or overlap fermions, incurring in 
the cost another factor of ~ lOOx; the correct light quark masses are ~ 2.5 MeV, 
and 5 MeV, not 20 MeV; the box size should be ample enough to accommodate 
the hadrons of interest (the Compton wavelength of the pion is about 1.4 fm, while the 
scattering length for the deuteron is about 5 fm); and the lattice spacing of the real world 



is, of course, zero. Finally, the above cost estimate only covers generation of lattice 
configurations; one must also account for the cost of generating quark propagators, the 
number of which grows factorially with the number of quarks involved — an unfortunate 
fact highly relevant to the study of even the smallest nuclei! Lattice QCD studies of a 
helium nucleus, for example, require 6!^ = 518,400 quark propagator contractions. It 
is easy to see that a brute force approach on a Pflops machine will not provide useful 
information about the a particle at realistic quark masses. 

Effective field theory is the tool that will allow us to extract useful information from 
available technology, giving us the ability to simulate real systems at unrealistic lattice 
parameters. 

In particular, we will have the opportunity to learn about fundamental properties of 
matter which are not directly obtainable from experiment, and which are necessary 
inputs for reliable nuclear structure or equation-of-state calculations. These include an 
improved understanding of three-body forces, such as in the experimentally inaccessible 
7 = 3/2 channel, and the interactions between hyperons and nucleons. The thesis of this 
talk is that progress in these directions will need an intense effort by theorists in order to 
extract physically relevant quantities from feasible lattice calculations, and that the basic 
tool for this effort will be effective field theory. 

Effective field theory (EFT) in all its forms is basically a perturbative expansion in 
the ratio of two length scales. As such, its validity requires that there be small ratios to 
exploit. Chiral perturbation theory has been the basic EFT exploited in continuum QCD, 
making use of the mass gap between the pion and the heavier hadrons. In addition, 
an effective field theory for nuclear physics has been in the making over the past 15 
years, which incorporates an additional small ratio, the QCD length divided by the A^A^ 
scattering length. 

What is new when working with lattice QCD is that there are a host of additional 
dimensionful scales which do not exist in the real world, but which can be profitably 
exploited. These include the lattice spacing, the lattice size, and independently varied 
masses for valence and sea quarks. EFT allows one to 

• extrapolate to smaller quark mass than is feasible to simulate 

• parametrize and correct for finite lattice spacing errors 

• parametrize and correct for finite volume errors 

• extract physics from "cheaper" fermions 

• determine 5-matrix elements from Euclidean simulations by measuring volume 
dependence of the spectrum 

• extract useful physical quantities from complicated multi-hadron systems 

USES OF CHIRAL PERTURBATION THEORY 
Quark mass extrapolation 

Chiral perturbation theory is an expansion of the Lagrangian for low energy QCD 

about the chiral limit, niq = 0. As such, it is obviously useful to extrapolate from lattice 
simulations at somewhat heavy quark masses, down to realistic quark masses. For this 



to work, the lattice quark mass has to be light enough so that the chiral expansion 
still converges. The chiral expansion parameter for mesonic processes is m|/A^, where 
A ~ Mp is not far from 1 GeV. A light quark mass of 20 MeV, for example, corresponding 
to rrijj: ~ 325 MeV, should be within the range of validity of chiral perturbation theory. 
I will not dwell on this conventional and important application of chiral perturbation 
theory which is widely familiar (see lectures by S. Sharpe [5] for a comprehensive 
introduction to lattice applications of chiral perturbation theory). 

Lattice spacing extrapolation 

Another application of chiral perturbation theory is to account for finite lattice spacing 
errors. One first matches the lattice action to the "Symanzik action" - a continuum 
theory with all operators allowed by the lattice symmetries, suppressed by powers of 
the lattice spacing a appropriate to the dimension of the operator. For example, with 
Wilson fermions (which do not possess a chiral symmetry), the leading operators in the 
Symanzik action not present in continuum QCD include 

• dimension-3 chiral symmetry violation: a^qq 

• dimension-5 chiral symmetry violation: aqo^yGfivl 

• dimension-6 Lorentz violation: a^qD^^j^q 

In order to determine the effects of finite lattice spacing on low energy hadronic 
physics, on can then match the Symanzik action onto a generalized chiral Lagrangian, 
which includes the effects of these finite lattice spacing operators [dlv|, S]. The coeffi- 
cients of these operators may be determined by making measurements at several lattice 
spacings, and then the extrapolation to a = may be improved. 

This program is versatile and can be applied to different lattice fermion formulations. 
For Wilson fermion the chiral symmetry violating operators give rise to the leading 
0{a) corrections, even when the dimension-3 operator is fine-tuned away. For staggered 
fermions, corrections begin at 0{a^), but the effective theory is complicated by the 
presence of additional "tastes", with an approximate SU (4) taste symmetry, broken by 
finite lattice spacing operators. The analysis of the chiral Lagrangian is simplest for 
chiral lattice fermions, such as domain wall or overlap fermions, which automatically 
avoid the 0{a) operators without incurring spurious fermion tastes. The computational 
price of dynamical chiral fermions is about a factor of 100, which is severe. 

Partially quenched chiral perturbation theory 

Quark masses appear in two distinct ways in the calculation of a correlation function 
in lattice QCD: either in the fermion determinant, which controls the gauge field con- 
figuration one generates; or in the fermion propagators on sews together in the gauge 
field background to compute the desired Green function. The former is called the "sea 
quark mass", the latter the "valence quark mass". In the real world they are the same, 
but in a lattice calculation they can be different. By making the valence quark mass light 




FIGURE 1. Schematic representation of parameter space in partially quenched theories, from ifTTll . 



while keeping the sea quark mass heavier, one obtains the benefits of chiral symmetry to 
leading order in the gauge coupling, without paying the light fermion price in generating 
the gauge field configurations. 

One can think of this unphysical partition function arising from unphysically heavy u 
and d "sea" quarks with mass nis, plus two flavors of light "valence" quarks with mass 
niy, plus two flavors of ghosts with mass my to cancel the valence quark contribution 
to the fermion determinant. Physical hadrons are then those made of VV quarks (as 
well as the strange quark), in the limit my = nis- The theory has additional unphysical 
mixed states of VS and SS content. A chiral Lagrangian can be constructed for this 
system, which contains new operators, but also all the operators of real QCD ^ [lo|] . 
Provided that ms is light enough for chiral perturbation theory to apply, by varying nis 
and my independently, and by making use of generalized chiral perturbation theory for 
the various types of correlation functions one can compute, it is possible to isolate and 
measure quantities of interest in QCD, such as the Gasser-Leutwyler coefficients for 
NLO chiral perturbation theory [11]. See fig. [T] for a picture of the expanded parameter 
space. For an example of how to use the partially quenched method to determine the 
up quark mass, see [fl2ll : there it was shown that by computing meson masses in the 
combination {Myy +^55 — 2M|y) one can extract the combination of Gasser-Leutwyler 
coefficients {2L% —^5) for QCD (where S and V label the two propagators used). 



Mixed action 

There has been much work done recently with staggered fermions, employing "the 
fourth root trick" to reduce the number of tastes from four to one per physical quark 
flavor. This results in a nonlocal theory at finite lattice spacing, and there has been a 
controversy about whether or not the resulting lattice theory is in the same universality 
class as QCD and is capable of delivering an approximation to continuum QCD [13l[l4ll. 
In addition, the EFT for staggered fermions is extremely complex and difficult to work 
with for baryons, due to the multiplicity of tastes and taste symmetry violating operators. 



The benefits of staggered fermions are their computational cheapness, and so at least for 
now, they are widely employed for mesons (see, for example, [,15,1 ). 

It is possible to improve the utility of staggered fermions by working with mixed 
actions, where the sea quarks are staggered, while the valence quarks are domain wall 
fermions. This approach benefits from combining the speed of staggered fermions and 
the chiral symmetry of domain wall fermions, and has been used in recent calculations 
of gA at the ~ 10% level [16], the ratio fn/fK to ~ 1% [[13], n-n scattering [lS],K-n 
scattering [fl9ll . and two-nucleon properties [|20ll . For mixed action computations one 
can use the partially quenched chiral perturbation machinery described above to extract 
physical results from the lattice calculations [21]. 



Volume dependence 

An additional handle on QCD provided by the lattice is the ability to manipulate the 
volume. Chiral perturbation theory is useful for understanding volume dependence of 
physical quantities, because it is the lightest modes that are most sensitive to finite vol- 
ume. For a sufficiently small lattice, the zero-momentum pion modes become collective 
coordinates corresponding to global rotations of the chiral condensate. This occurs when 
m{qq)V < 1, defining the e-regime [22]. This regime was recently cleverly exploited to 
extract information about QCD in the infinite volume continuum for the A/ = 1/2 rule 
[|23|] . and for nucleon properties [ | 24|] . Both references make extensive use of chiral per- 
turbation theory. 



NUCLEAR EFFECTIVE FIELD THEORY 

EFT will play a critical role in computing nuclear physics properties on the lattice. Just 
as for meson interactions, one would like to find strategies to measure on the lattice the 
coefficients of the most relevant operators which control the interactions of nucleons, 
and then use that effective theory to compute nuclear properties. This program requires 
(i) that there be a sensible EFT for the interactions of nucleons, and (ii) that one can 
relate lattice measurements in Euclidean space to experimentally measurable quantities, 
such as scattering lengths. 



Nuclear effective theory was pioneered by Weinberg [|25|| . first explored in ref. [1260, 



and then developed by many subsequent authors (for a review, see 112711 '). For low en- 
ergy nucleon interactions, t he p ion may be considered as heavy, and the pion-less 
EFT developed in S S, Efl '3l'] consists of contact interactions ~ Co{N'^Nf + 
C2iN'^N){N'^V^N) + 0{p^) . . . The momentum expansion treats C2„ = when 
renormalized at a scale = 0{p), and Q is summed to all orders as the leading con- 
tribution, while the higher €2,1 are inserted perturbatively. For p > niTi/l, the pion has 
to be included explicitly. Unfortunately, unlike the case for chiral perturbation theory, 
the power counting scheme for the effective theory for nucleons interacting via pions is 
somewhat controversial. The reason is in part because NN scattering is nonperturbative, 
and so the actual scaling of an operator does not match its naive dimension, making it 
difficult to construct a consistent power counting scheme. The original Weinberg scheme 



suffers inconsistencies, where counterterms appear at higher orders than the divergences 
they are supposed to cancel. One example of this was given in [32] where it was shown 
that a quark mass-dependent counterterm for the non-derivative A^A^ vertex was required 
at leading order; a more recent numerical analysis of NN scattering by the tensor force 
demonstrated that counterterms are needed at leading order in all partial waves where 
the interaction is attractive [33], even though such counterterms would be subleading in 
Weinberg's expansion. By working at a fixed and not too large cutoff, this problem can 
be swept under the rug, but this procedure in effect corresponds to constructing a model 
for short distance physics rather than performing a bona fide EFT calculation. This point 
of view is not universally accepted, and for recent contributions on various sides of the 
controversy see refs. [34, 35, 36]. The KSW expansion [29, 30] was offered as an alter- 
native to Weinberg's expansion, but was found not to converge well for two nucleons in 
the = 1 channel [|37l] . I believe that the most consistent expansion currently available 
is that of ref. fs^, with generalizations to account for three nucleon forces [39]. 

There have been many notable successes of the nuclear effective theory, most remark- 
ably at very low energy. A nice example is the isolation of the EFT coupling Lia, the 
analogue of gA for the axial isovector two-nucleon current. By fitting it to data one can 
compute the neutrino-deuteron breakup cross section, thereby reducing a major source 
of systematic error in the analysis of data from the Sudbury Neutrino Observatory po( ]. 

In order to study multi-nucleon states on the lattice, one approach is to compare EFT 
and numerical results in Euclidean space, and determine the EFT couplings. Another 
approach is to compare lattice results for S -matrix elements with the predictions of 
the EFT. An important contribution to the problem of extracting S-matrix elements 
from Euclidean lattice theory was devised by Liischer, who showed how the volume 
dependence of the energy for a 2-particle state in a box yields the scattering lengths 
[41|. A nonrelativistic formulation found in ref. |4^|2^ goes as follows: The Feynman 



scattering amplitude for two nucleons has the form £/ = {An/ M) / {pcoi5{p) —ip), 
where p = \JME, E being the energy in the center of mass. When formulated in a box, 
the energy eigenvalues correspond to zeros of Re[{sz/)^'^]. So the energy eigenvalues 
solve 

= i?e[K)-i]box=^e[K)-i]L=o,+ (i?e[K)-^]box-^e[(^/)-i]L=oc) . (2) 

The first term on the right is just proportional to pcot5 in an infinite box. The second 
term is the difference between the bubble diagram for two nucleons scattering off each 
other computed in an infinite box versus computed in a finite box. This is a finite and 
computable function of the box size L and the eigenenergy £„, and one arrives at the 
formula 

1 ^ , 

Pn cot 8 (pn) = —S{ri„) , p^ = VEjl , 7]„ = {pnL/lTtY , (3) 



kL 



5(77) = lim 



(4) 



where the j are integer triplets. By measuring the energy eigenvalues for two nucleons 
in a box, one can then in principle solve the above equation for 5{p„). A pioneering 



measurement of two nucleon scattering lengths using this method was performed in 
ref. [|201 . where they concluded that one would need a box of size 5 — 15 fm in order to 
study properties of the deuteron. This is a large box, but much smaller than one would 
conclude directly for Liischer's work, which would seem to indicate a. 

It is important to push two nucleon studies much further, technically and theoretically. 
Using EFT techniques, one should learn how to measure matrix elements of currents in 
the two-nucleon state, study hyperon-nucleon interactions, and prepare the groundwork 
for the study of three-nucleon states on Pflops machines. 

THE FRONTIER 

Lattice QCD can make substantial and useful contributions by predicting properties of 
baryons which are not experimentally accessible. I believe that in the foreseeable future, 
useful prediction are feasible in four areas: 

• The masses and couplings of QCD resonances and hybrids; 

• Strangeness physics, such as ^A'^, YY and YN interactions, where K is the anti-kaon, 
y is a hyperon, and is a nucleon; 

• Determination of 3-body interactions, such as in the 7 = 3/2 channel; 

• Quark mass dependence of nuclear properties. 

The first category could make an important contribution to the JLab experimental pro- 
gram; information the second category could answer basic questions about dense matter 
in neutron stars, such as which hadronic channel is favored when strange quarks first ap- 
pear, the K , A, or E ? Accomplishing the third task would be an important milestone, 
whereby lattice QCD could inform the so-called ab initio nuclear structure calculations 
and make them significantly more ab initio. The fourth category could be of interest in 
understanding how fine-tuned is our world, and could be important in certain cosmolog- 
ical theories where the quark masses are dynamically determined quantities [j43l] . 

All of these projects could be very rewarding, and will require an intensive theoretical 
effort that further develops the 3-nucleon EFT, extends Liischer's work to three particles 
in a box (where inelastic thresholds could cause problems 10] and improves the 
available computational algorithms. Lattice QCD will clearly play an important role in 
the future progress of nuclear theory, and EFT will be a vital component of the program. 
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